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ABSTRACT 

Context. Thanks to recent large scale surveys in the near infrared such as 2MASS, the galactic plane that most suffers from extinction 
is revealed and its overall structure can be studied. 

Aims. This work aims at constraining the structure of the Milky Way external disc as seen in 2MASS data, and in particular the warp. 
Methods. We use the Two Micron All Sky Survey (hereafter 2MASS) along with the Stellar Population Synthesis Model of the 
Galaxy, developed in Besancon, to constrain the external disc parameters such as its scale length, its cutoff radius, and the slope of 
the warp. In order to properly interpret the observations, the simulated stars are reddened using a three dimensional extinction map. 
The shape of the stellar warp is then compared with previous results and with similar structures in gas and dust. 
Results. We find new constraints on the stellar disc, which is shown to be asymmetrical, similar to observations of HI. The positive 
longitude side is found to be easily modelled with a S shape warp but with a slope significantly smaller than the slope seen in the HI 
warp. At negative longitudes, the disc presents peculiarities which are not well reproduced by any simple model. Finally, comparing 
with the warp seen in the dust, it seems to follow a slope intermediate between the gas and the stars. 
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1. Introduction 

In recent years wide surveys have been obtained thanks to wide- 
field mosaic of CCD cameras and helped by the availability 
of dedicated telescopes both from the ground and from space. 
These large data sets benefit Galactic structure studies. 

However, interstellar extinction remains a serious obstacle 
for the observation of stars in the Milky Way, and for interpret- 
ing these observations in terms of Galactic structure. As the ex- 
tinction suffered in the near infrared is roughly an order of mag- 
nitude lower than in the visible, recent infrared surveys provide a 
suitable tool to probe further into the densest parts of the Galaxy, 
namely the Galactic plane. 



Among them, the 2MASS survey ( [Skrutskie et al.[|2006| ) has 
provided one of these rich data sets from which substantial re- 
sults have been obtained. The homogeneity of these data sets 
is a great help for avoiding systematic bias and bad links from 
one set to another, which have created troubles in data interpre- 
tation in the past, disentangling real cosmic substructure from 
local zero point errors for example. 

The analysis of near infrared data in the whole Galactic plane 
can bring important clues regarding the stellar populations. In 
particular, a good knowledge of the disc structure is crucial to 
characterize the overdensities found close to the Galactic plane, 
such as the Canis Major overdensity. This structure locate d at 
/ = 240°, b = -8° was revealed by |Martin efaT] ( |2004a| ). Its 
nature is still under debate. While several studies based on kine- 
matics, star count s and CMP studies of this field conclude that it 
is a dwarf galaxy ([Martin et al.[[2004bl|2005t[Martinez-Delgado| 
eTaI^[2005l|Butler||2007t|De Jong| |2007| ), other authors tested 
several warp models and claimed that the warp is still a possi- 
ble explanation of the Canis Major overdensity (Momany et al. 



2004] |2006[ |L6pez-Corredoira| 12006) . It shows how important it 
is to have a good description of the Galactic disc structure. 



Warps are a common feature in external galaxies (Bosma 



1991). Radio observations revealed that the Milky Way itself 
hosts a warp and fl ared disc. The warp and flares have been ob- 
served in HI data (Hender son et al.| |1982| |Burton & te Lintel 
Hekkert||1986l|Burton||1988t Diplas & Savage||199T 



& Sofue[ 120031 JLevine et al.| 12006]), in mole cu 



Nakanishi 



ar clouds 



(Grabelsky ~et al.| |1987[ |Wouterlout etal) [1990} |May et al. 



|19971 ), from OB stars ( [Miyamoto et al.|| 1988 ||Reed[Tl 996) and 



final ly in CQBE data jPorcel et al.| |1997| |Freudenreich| |1 998 
Drimmel & Spergel, 2001). The warp and flare in the Galaxy 



ha ve also been detected in the i n terstellar dus t (Freudenrei ch et 
[aTl [T994l|Marshall et al.|[2006]).|Derriere & Robin| 42001) ) have 
analysed DENIS (Epchtei n et al.||1999| ) star counts in several 
fields in the Galactic plane and showed evidence of a warp and 
a flare in the stellar disc. More recently Lopez-Corredoira et al. 
(2002) modelled a stellar warp from the analysis of the stellar 
populations in the 2MASS survey. 

Warps may have originated from interactions between the 
disc and (i) the dark halo (if angular momenta are not aligned), 
(ii) nearby satellite galaxies, such as the Sagittarius dwarf or 
the Magellanic clouds, (iii) infalli ng intergalacti c gas. From the 
analysis of their angular momenta, Bailin ( 2003 ) argues that the 
Milky Way warp and the Sagittarius Dwarf Galaxy may be cou- 
pled. |Garcia-Ruiz (2002) used N-body simulations to explore 
the possibility that the warp has been created by tidal forces 
caused by the Magellanic Clouds. They find that neither orienta- 
tion nor ampl itude of the warp can be reproduced this way. On 
the contrary, [Weinberg & Blitz ( 2006 ) using perturbation the- 
ory showed that a Magellanic Cloud origin for the warp of the 
Milky Way can explain most quantitative features of the outer HI 
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layer. Alternatively Lopez-Corredoira et al.| ( 2002| ) propose that 
the warp is due to intergalactic accretion flows onto the Milky 
Way disc. 

Here we report on an the analysis of the external stellar disc 
of the Galaxy, and particularly the warp, as seen from near in- 
frared star counts using a Galaxy model. We overview the stellar 
population synthesis approach and the basic inputs of the model, 
in particular for the external disc (Sect. [2]). In Sect. [3] we briefly 
describe the determination of a 3D map of interstellar extinction 
and show a first comparison between the Galaxy model and the 
2MASS survey in the Galactic plane. We then describe (sect. [4]) 
constraints obtained from the comparisons of model predictions 
using different disc models with 2MASS data and compare this 
with results obtained from other studies. In sect. [5] we provision- 
ally conclude on the comparison between the gaseous warp and 
the stellar warp and about the origin of this warp. 



2. The Galaxy model 

2.1. Stellar population synthesis approach 

In order to analyse the large data set provided by 2MASS, we 
use a stellar popul ation synthesis approach (the Galaxy model, 
Rob in et al. ( 2003| )). This model aims at assembling current sce- 
narii of galaxy formation and evolution, theories of stellar for- 
mation and evolution, models of stellar atmospheres and dynam- 
ical constraints, in order to make a consistent picture explaining 
currently available observations. When scenarii of Galaxy for- 
mation and evolution are inferred from suitable constraints, they 
may be tested using population synthesis models whose predic- 
tions can be directly compared with observations. This synthetic 
approach ensures that biases have been correctly taken into ac- 
count and that the scenario is compatible with many kinds of 
constraint. 

The validity of any Galactic model is always questionable, as 
it describes a smooth Galaxy, while inhomogeneities exist, either 
in the disc or the halo. The issue is not to make a perfect model 
that reproduces the known Galaxy at any scale. Rather one aims 
to produce a useful tool to compute the probable stellar content 
of large data sets and therefore to test the usefulness of such 
data to answer a given question relative to Galactic structure and 
evolution. 

The model we use is the Galaxy model developed in 
Besancon. Its originality, as compared to a few other popula- 
tion synthesis models presently available for the Galaxy, is the 
dynamical self-consistency. The Boltzmann equation allows the 
scale height of an isothermal and relaxed population to be con- 
strained by its velocity dispersion and the Galactic potential 
( [Bienayme et al. 1987| ). The use of this dynamical constraint 
avoids a set of free parameters quite difficult to determine: the 
scale height of the thin disc at different ages. 



2.2. Basic scheme 

The model assumes that stars belong to four main populations: 
the thin disc, the thick disc, the stellar halo (or spheroid), and 
the outer bulge. The modelling of each population is based on a 
set of evolutionary tracks, assumptions on density distributions, 
constrained either by dynamical considerations or by empirical 
data, and guided by a scenario of formation and evolution, that 
is to say assumptions on the initial mass function (IMF) and the 
star formation rate (SFR) for each population. More detailed de- 
scriptions on these con straints can be found i n |Haywood et al.| 
(TWJ\ for the thin disc, |Reyle & Robin| ( [2QQT] ) for the thick disc, 



Robin et al.| ( [2Q00] ) for the spheroid, and |Picaud & Robin| ( [2004| ) 
for the outer bulge. Recently, the Hipparcos mission and large 
scale surveys in the optical and the near-infrared have led to 
new physical constraints improving our knowledge of the over- 
all structure and evolution of the Galaxy. These new constraints 
are included in the version of the model described in lRobin et al.l 
( [2003] ). 

2.3. The external disc model 



The disc is assumed to be axisymmetric and spiral arms are 
not modelled. The thin disc is divided into 7 age components. 
The distribution of each disc component is described by an 
axisymmetric ellipsoid with an axis ratio depending on the 
age, using an Einasto law. The disc scale length, h R =2530 
pc, is derived from the fitting of the old disc parameters on 
DENIS star coun ts towards the Galactic bulge direction ( |Picaud| 
& Robin |2004| ). Several studies have shown that the edge of 
the disc is detected at a galactocentric distance of about 14 kpc 
( [Robin et aL] [T992 ; Rup hy et aL| [1996 ). The disc is also warped 
and flared and all stellar components in the disc have the same 
warp and flare. 

As in |Gyuk et al. ( |1999| ), we model the flaring by increasing 
the scale heights by a factor &fl are , beyond a galactocentric radius 

^flare- 

^flareCR) = 1 + 7flare(^ _ ^flare) 



The amplitude is taken from Gyuk et al. ( 1999) (yfl are = 5.4 x 
10" 5 pc" 1 ) and the radius 7?fl are at which the disc starts flar- 
ing from Derriere & Robin ( [2001 ).They found using DENIS 
data that the minimum radius of the flare might depend on 
the longitude considered and determined an average value of 

tfflare = 9.5 kpc. 

We first model the warp as a S shape warp, with a tilted ring 



model following Porcel et al. (1997]). The height z warp of the warp 



over the plane b=0 u as a function of galactocentric coordinates 
is computed using: 



ZwarpW = Twarp X (R - # warp ) X shl(0 - 



Pwarp 



.) 



Following Burton (1988), we assume, as a starting point, that 
the Sun lies on the line of nodes of the warp (0 warp = 0°). Most 
of warp studies indicate values of this angle close to this one 
within an uncertainty of a few degrees. We use R warp = 8.4 kpc 
indicated as th e best value for the startin g galactocentric radius 
of the warp by Derr iere & Robin| ( |2001| ). This value is close to 
the galactocentric distance of the Sun as adopted in the Galaxy 
model. We adopt, as an initial guess, the value of the displace- 
ment of the mid plane y warp = 0.18 from Gyuk e t"aD(|1999|). It 
gives positions of the mid-plane similar to the Lopez-Corredoira 



|et al.| ( [20Q2] ) andjprimmel & Spergel (2001) stellar warp models, 
with a maximum height of 630 pc at 12 kpc from the Galactic 
Centre. 



3. Interstellar extinction 

The 2MASS survey allows us to study large scale structure in the 
Galaxy by comparing with simulations of the Galaxy model, par- 
ticularly in the Galactic plane where extinction is quite high. The 
comparison is realistic only if the extinction is known. Without 
a good estimate of the extinction and of its distance it is nearly 
impossible to understand the structure in the thin disc. Hence it 
is imperative to start by constructing a three dimensional extinc- 
tion map of the Galaxy. 
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Extinction is so clumpy in the Galactic plane that it de- 
termines for a great part the number density of stars, more 
than any other large scale stellar structure. Thus, it is possi- 
ble to extract information about the dis tribution of the extinc- 
tion from photometry and star counts. Marsh all et al.| ( 2006| ) 



have shown that the 3D extinction distribution can be inferred 
from the colour distributions in the 2MASS survey. Using stellar 
colours in J - Ks as extinction indicators and assuming that all 
of the model prediction deviations on small scales from observed 
colours arises from the variation of extinction along the line of 
sight, they built a 3D extinction map of the galactic plane. The 
final resolution in longitude and latitude is 15 arcmin and the 
resolution in distance varies between 100 pc to 1 kpc, depend- 
ing on stellar density and on the dust distribution along the line 
of sight. The resulting 3D extinction map provides an accurate 
description of the large scale structure of the disc of dust. These 
maps clearly show the dust warp in the external disc (Fig. 1). 
Note that the warp found for the dust differs from the stellar 
warp model. 

After applying the extinction to the simulated stars, it is pos- 
sible to perform a realistic comparison of the number n mo ^ of 
stars from the Galactic model with the number n \> s of observed 
stars in the 2MASS data (Fig. [2]). The comparison points to sig- 
nificant discrepancies. Thus it appears that several galactic com- 
ponents are not properly characterized in the Galaxy model: (i) 
the bulge extends too far from the plane, (ii) the disc in the inner 
part of the Galaxy (|1| <75°) has a scale height slightly too large, 
(iii) the warp of the external disc is too strong, and (iv) the mod- 
elled warp is symmetrical whereas the observations point to an 
asymmetrical warp. 

In this work we concentrate on the external disc (|1| >60°) 
and propose to derive a better parametrization in three steps: (i) 
simulations of the Galaxy are obtained using different parame- 
ters, (ii) as the extinction found depends on the parameters 
in the Galaxy model (the sensitivity of the extinction mod el 
to the Galaxy model is discussed in Mars halTet al.| ( [20 06 }), 
the extinction is derived for each set of Galactic parameters 
using the Mars hall et al.| ( |2006| ) method and applied to the sim- 
ulated stars, (iii) number counts of simulated and observed stars 
are statistically compared. 

4. Results 

As seen in Fig. [2| the assumed slope of the warp is not suitable 
for correctly reproducing the data. The slope value is too high 
and the simulated warp deviates too far from the b = 0° plane. 
Too many stars are predicted below the galactic plane at nega- 
tive longitudes and above the plane at positive longitudes. The 
overestimate is higher for / < 0°, showing that the stellar warp 
is asymmetric. This can be either an apparent asymmetry if 
the Sun does not lie on the line of nodes, or a real asymme- 
try in the disc. However the effect of the position of line of 
nodes does not solve the problem, as explained further be- 
low. Another possibility is that the asymmetry is due to the 
distribution of the dust with respect to the stars. In our ap- 
proach we assumed as a first approximation that a symmet- 
ric warp exists in the stars only. The extinction we derive 
does not suppose any geometry for the large scale distribu- 
tion of the dust - it simply finds the extinction necessary to 
correct the simulated stellar colours to fit the observed ones. 
The dust warp then appears naturally in the extinction dis- 
tribution and this dust warp is not symmetrical. Assuming 
this dust warp, we then tentatively ajust the star counts on 
both sides of the galaxy, in order to check the values of the 



northern and the southern warps independently. As such, it 
is not possible for us to create a stellar overdensity in the 
model by a misalignment of the stellar and dusty warps. 

Fig. [3] (top panel) shows 2MASS star counts in the longi- 
tude ranges 120° to 60° and -60° to -120°. They are compared 
with simulated star counts obtained with the simulations from 
the Galaxy model (middle panel). The difference is shown in the 
bottom panel. In addition to the overestimate of the warp slope, it 
appears that the modelled density of stars is too large at |/| > 80° 
probably due to the value of the scale length h R being too high. 

A new simulation with the smaller values y W arp — 

0.09 and 

Hr = 2200 pc is shown in Fig [4] It substantially improves the 
agreement between model and data at positive longitudes. The 
value h R =2200 p c is in agreement with the one determined by 
|Robin et aL| ([T992) from star counts towards the anticenter. 

However the comparison is still unsatisfactory at negative 
longitudes. There remains a significant discrepancy, which can 
not be improved by varying the assumed slope, the disc scale 
length or the sun-centre distance (assumed to be here 7.9 kpc). 
Further testing by varying the line of nodes parameter warp 
(from -30° to 30°) shows that it cannot account for this disc 
asymmetry either (for the extreme values, the fit is even de- 
graded). 

At negative longitude, this model parametrization (S shapoe 
warp) is clearly unsatisfactory. In an attempt to improve the fit 
at negative longitudes, we investigate other warp models, in par- 
ticular non-axisymmetric ones. It is well known that the HI warp 
is asymmetric: the gas warp bends back to the Galactic plane in 
the southern hemisphere at R>15 kpc (Burton 1988). The stellar 
warp may also follow this feature. 

Levine etaT] ( 2006 ) studied in detail the warp structure from 



HI surveys and found that the galactic HI warp is well described 
by a vertical offset plus two Fourier modes of frequency m-\ 
and 2, the m=2 mode accounting for the asymmetry. The height 
of the mid plane from the b = 0° plane takes the following form 
as a function of the azimuth cp : 

z(4>) = zo + z\sin((p - (pi) + Z2sin(2<p - <f> 2 ) 

where zo, Z\ and zi are the heights of the mid plane, quadratic 
functions of the galactocentric radius, <pi and 02 are the node 
lines of the Fourier modes 1 and 2, all of which are adjusted 
parameters. 

With the best fit parameters from Levine et al. (2006]) the 
m=l mode dominates up to about R<15 kpc. Closer to the 
Galactic centre the other modes are unimportant. The m=2 mode 
is significant at R>20 kpc. This means that the detection of these 
modes might not be easy with the 2MASS survey if the stellar 
warp follows the gas warp. Red clump giants at Ks = 14.3 (the 
limiting magnitude of 2MASS survey) reach distances of only 
about 10 kpc from the sun. 

In order to test whether the stellar populations might also fol- 
low this gaseous warp, we tentatively modelled the stellar warp 
with the |Levine et al.| ([2006) model. If we adopt their parame- 
ters, the simulated counts do not give an acceptable fit. However 
if we adopt a closer value for the starting galactocentric radius of 
the m=2 mode, 12 kpc instead of 15 kpc, and reduced slopes, the 
fit at negative longitude is slightly better and still good at posi- 
tive longitude, although not as good as the tilted ring model with 
7warp=0.09 (Fig.|4j). The result is shown in Fig. 5. The systematic 
discrepancy remains which is not explained by the two Fourier 
modes warp model. 

We also investigated the possibility of a varying value for 
the minimum of the flare, such as proposed by Derriere & Robin 
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Fig. 1. Warp in the dust as seen from 



Marshall et al. 



(2006) 3D extinction map. The plot shows the density of absorbing material 
at 8kpc from the sun, as a function of longitude and latitude. The dashed line is the position of the mid-plane expected from a 
0.18 slope warp model, and the solid line is the position from the fitted slope of 0.14 at positive longitudes and 0.11 at negative 
longitudes. 



( 2001| ). However, a lower value 7?fl are = 8.5 kpc does not improve 
the comparison. 

A possible explanation is that the truncation in the stellar disc 
occurs at shorter distances on this side of the disc, since the stan- 
dard truncation used on the other side of the disc (R c = 14 kpc) 
creates an overestimation of the star counts at large distances. 
Indeed, considering a cutoff radius R c = 12 kpc slightly reduces 
the discrepancies at negative longitudes but still not satisfying 
(see Fig. 6). Assuming a lower value for R c causes a deficit of 
model stars around longitudes / = 70°. At positive longitudes, 
the model with a disc truncation at 12 kpc does not agree as well 
as a model with R c = 14 kpc. This could be an indication of a 
varying distance of the edge of the disc as a function of the az- 
imuth. The edge of the disc possibly follows an ellipsoidal shape 
instead of a circular one. 

The mid plane of stars as a function of longitude is given by: 



Hi) 



Z b b.p(Lb) 



where p(/, b) is the density of stars. Fig. [7] shows the mid plane 
for the data (solid thin line) and the simulations (solid thick line) 
obtained with y warp = 0.09, Hr = 2200 pc. As expected from 
Fig. [4] the mid plane of simulated stars is in better agreement 
with that of data at positive longitudes. However, the mid plane 
of simulated stars is not able to fit the observed one over the 
whole longitude range. We also computed the mid plane of giant 
stars only (dotted line). The selection of giants in 2MASS data 
is based on their expected colour given the derived extinc- 
tion. The giants being the brightest stars probe a more distant 
region of the disc. For that reason, the mid plane of giant stars 
is further away from the b = 0° plane. At negative longitudes, 
the slope y warp = 0.09 results in a good fit of the mid plane of 
giant stars. However, at positive longitudes, this slope seems to 
be underestimated when considering giant stars only. This could 
be an indication of a varying slope of the warp as a function 
of galactocentric distance. Such a model with z warp vary ing with 
galactocentr ic distance has been proposed by |L6pez^C orredoira 
|et al.| ( |2002| ), based on the analysis of 2MASS data. However, 
the range of distance over which the warp is constrained is rather 
small, thus the observed variation in slope is hardly significant. 
Another interpretation can be that the warp depends on the 
age. If young stars are born in a gaseous warp originating 
from gas dynamics, only young stars will follow the warp. 
However, if the warp is the response to a perturbation by a 



dwarf galaxy or a misalignement of the dark halo, all com- 
ponents should be perturbed independently of their age. 



5. Discussion 

The warp is a structure that appears in all components of the 
Milky Way disc: stellar, dust and gas. Hereafter we review recent 
studies of these different warps and compare them. 



5.1. The stellar warp 

The maximum height of the mid plane at a galactocentric dis- 
tance of 12 kpc is 315 pc and smaller than the one found in other 
studies. 

The three-dimensional model for the Milky Way fitted to the 
far-infrared and near-infrared data from the COBE/DIRBE in- 
strument by Drimm el & Spe rgel ( 2001) assumes that the dis- 
placement of the warp follows a quadratic function of the galac- 
tocentric radius, z war p(R) = 27.4(7? - R warp ) 2 sin(0), with R warp 
= 1 kpc. The maximum height of the plane at 12 kpc from the 
Galactic center is 685 pc. 

With an elevation of the stellar warp varying as a power law 
as a function of the gala ctocentric distance, z wdLW (R) = 1.2 x 
10" 3 #(kpc) 5 25 sin(0 + 5°), |L6pez-Corredoira et at| ( |2002| ) found 



a similar value (556 pc), based on the analysis of the old stellar 
population using data from the 2MASS survey. 

Momany et al.| ( [20 06 ) derived the structure of the stellar warp 



using 2MASS red clump and red giant stars selected at helio- 
centric distances of 3, 7 and 17 kpc. Their results are in agree- 



ment with Drimmel & Spergel (2001 ) and Lopez-Corredoira et 
al. (2002]) results for the Southern warp, whereas the Northern 
warp (at positive longitudes ) is less p ronounced. Their results are 
well in agreement with the Yusifov (2004 ) model based on the 
analysis of the distribution of pulsars. They invoked two factors 
that may contribute to the appearance of an as ymmetric warp . 
First the Sun is not located on the line of nodes ( Yusifov (2004) 
found a rather high value warp = 15°). Second the Northern 
warp is located just behind the Norma-Cygnus arm and possible 
variations of the extinction may produce an apparent asymmetric 
warp. 

Fig. 8 shows the height of the galactic disc as a function of 
longitude at an heliocentric distance of 7 kpc (comparable to 
Fig. 11 in |Momany et al.| ( [2006| )) from these studies of the stel- 
lar warp. On the contrary, our model of the Northern warp is in 
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agreement with Moman y et al. (2006 ) but our Southern warp is 
smaller than in all other studies. 




300 200 100 

Galactic longitude (°) 



Fig. 8. Elevation of the disc mid plane as a function of longi 
tude models at an heliocentric distance of 7 kpc from different 
studies. Solid line: this work y warp = 0.09, dashed line: |L6pez- 
Corredoira et al. 1(2002), dotted line: |Drimmel & SpergeipnflT] ); 
squares: |Momany et al.| ( |2006] ). 



5.2. The dust warp 



Marshall et al.| ([2006 ) detected the dust warp when convert- 
ing the 3D extinction map into interstellar dust volume density. 
Fig. 1 shows the distribution of dust at 8 kpc from the Sun as 
a function of longitude. The warp is clearly seen on both sides, 
but with significant differences. They determined the slope of 
this dust warp to be 0.14 at positive longitudes and 0.11 at neg- 
ative ones. They also noticed that the warp seems to start earlier 
at negative longitudes (at 7.8 kpc rather than 8.9 kpc found on 
the other side). 

The displacement at a galactocentric distance of 12 kpc is 
460 pc, larger than the 315 pc value found for the stellar warp in 
this work. 



5.3. The gas warp 



Nakanishi & Sofue (2003 ) found that the HI disc is warped and 
asymmetric, with a displacement larger a t positive longitudes, in 
agreement with the dust warp found by Mars halTet al.| ( |2006| ) 
and the stellar warp in this work. However, the HI warp starts at 
a galactocentric distance of 12 kpc, which is close to the edge of 
the stellar disc. 

Recently, [Levine et aT| ( |2006| ) used a superposition of Fourier 
modes to describe the HI disc. The m= 1 mode is dominant at dis- 
tances smaller than 15 kpc. They found that this mode is nearly 
linear from 10 to 15 kpc, corresponding to a slope of 0.20, much 
higher than what is found in the dust and stellar warps. 
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Fig. 9. Maximum elevation of the disc mid plane from dif- 
ferent warp models. Thick solid line: HI |Levine et al.| (2006) 
m=l mode, thin solid line: stellar warp (this paper), dotted line: 
Marshall et al. ( 2 006]) dust model at positive longitudes, dashed 
line: Marshall et al.|(|2006) dust warp at negative longitudes. 



line of nodes. A comparison of the warp elevation as measured in 
dust, HI gas and stars is given in Fig. 9 together with the galacto- 
centric distance at which the warp starts. In our study, the stellar 
warp slope is significantly smaller than the HI gas warp, by a 
factor of about two. There is an indication that the HI gas is the 
strongest, the dust warp seems a bit smaller and the stellar warp 
is significantly smaller. However if we limit the investigation to 
short distances (R<12 kpc) the differences between various warp 
models are less important. It is probable that the difference we 
find with alternate models ([Drimmel & Spergel 200 1| Lopez- 
|Corredoira et al. 2002| Momany et al. 2006| ) originates from 
remote counts at galactocentric distances between 12 and 14 kpc. 
This comparison tends to show that the different components re- 
act differently to the forces at the origin of the warp. 

In external galaxies, |van d er Kruit ( 2008| ) noted that stel- 
lar discs look flatter than gas layers. This is understandable in a 
scheme where the HI warps start close to the truncation radius 
(truncation seen in the exponential distribution of stars which 
may be due to a threshold effect in the star formation efficiency). 
In this case the stellar warp is only visible in a small range of 
galactocentric radii. Robin et aT| ([l992) measured the truncation 
radius in our galaxy to be at about 14 kpc. This is a bit larger than 
where the HI warp starts, apparently at 10 to 12 kpc, and even 
larger than where the stellar warp starts (at ~ 8 kpc galactocen- 
tric radius). Hence the stellar warp is expected to be detected 
only between 8 and 14 kpc, at most. 

|van der Kruit (2008 ) suggests a scenario where the inner disc 
formed initially in a rigid and flat structure, while the HI warp 
resulted from later infall of gas with different orientations of an- 
gular momentum (in the case of a dark matter halo not aligned 
with the one of the inner Galaxy). This on-going infall could be 
identified with the HI flare. 



5.4. Comparison 

These studies reveal a warp present in all components of the disc 
(stars, dust, gas), all of which are asymmetric and with a similar 



6. Conclusion 

From the comparison of 2MASS star counts with a model of 
population synthesis in the outer regions of the Galaxy we in- 
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vestigate the external disc structure and in particular the warp 
feature followed by stars. We obtain a good fit of one side of the 
disc with a scale length ^=2200 pc and a simple warp model 
following a slope of 0.09 and starting at 8.4 kpc. On the other 
side, at negative longitudes, the fit remains unsatisfactory either 
with a S shape warp, or a Fourier modes shape. A shorter trun- 
cation of the disc improves a bit the fitting at negative longitudes 
but degrades it at positive longitudes. This could be an indication 
of a variation of the distance at which the stellar disc stops as a 
function of azimuth. Despite the uncertainties, it is clear from 
our study of the 2MASS star counts that the warp is less marked 
in stars than in the gas and it also shows different shapes of the 
disc at negative and positive longitudes. There is also a slight 
tendancy for the dust warp to be intermediate between the HI 
and the stars. Further studies are needed to confirm these con- 
clusions, in particular a dynamical approach with kinematical 
data, helpful in understanding the dynamical evolution of this 
structure, its origin and its link with the gaseous warp. 
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Fig. 2. Comparison between the number of stars from the Galaxy model n mo ^ and 2MASS data n ^ s . The counts are given in a colour 
table in number of stars per square degree to magnitude Ks = 12. The relative difference is given as (ft m od _ ^obs)/^obs- It varies 
from a factor of two excess in modelled star counts to a factor of two deficiency in modelled star counts. Regions where the Galaxy 
model overestimates or underestimates de density of stars clearly appear. 
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Fig. 3. 2MASS star counts (top), modelled ones (middle), and relative difference between the two (bottom) with y warp =0.18 and 
/z#=2530 pc. Dashed line indicates the b = 0° plane. On the left the Northern warp, on the right the Southern warp. 



8 



C. Reyle et al.: The Milky Way's external disc constrained by 2MASS star counts 





Fig. 4. Same as Fig. [3] with y warp =0.09 and /z#=2200 pc. The model which reasonnably fits the positive longitude side significantly 
differs from the data at the negative longitudes, showing a clear disymmetry between both sides of the Galactic disc. 





Fig. 5. Same as Fig. [3] with the two Fourier modes warp model from |Levine et al. ( 2006 ). This model does not better fit the stellar 
data than the simple model and the main discrepancy at negative longitudes is not removed. 
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Fig. 6. Same as Fig. |3] with a truncation of the disc at 12 kpc. The main discrepancy at negative longitudes is sightly reduced but not 
removed. 




Fig. 7. Mid plane of data (thin lines) and model (thick lines) for all stars (solid lines) and giants only (dotted lines). The simulations 
are obtained with y warp =0.09 and /z#=2200 pc. 



